Rabi-like oscillations of an anharmonic oscillator: classical versus quantum 

interpretation 
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We have observed Rabi-like oscillations in a current-biased dc SQUID presenting enhanced coher- 
ence times compared to our previous realization |py]. This Josephson device behaves as an anhar- 
monic oscillator which can be driven into a coherent superposition of quantum states by resonant 
microwave flux pulses. Increasing the microwave amplitude, we study the evolution of the Rabi 
frequency from the 2-level regime to the regime of multilevel dynamics. When up to 3 levels are 
involved, the Rabi frequency is a clear signature of quantum behavior. At higher excitation ampli- 
tude, classical and quantum predictions for the Rabi frequency converge. This result is discussed in 
the light of a calculation of the Wigner function. In particular, our analysis shows that pronounced 
quantum interferences always appear in the course of the Rabi-like oscillations. 

PACS numbers: ... 



I. INTRODUCTION 

It is well-known that the classical dynamics of an an- 
harmonic oscillator driven by an external monochromatic 
force can exhibit very complex and rich behavior, such 
as nonlinear resonances, bistability oscillations or chaotic 
dynamics [2, 3, 4, 5]. The search for a quantum signature 
in the quantum counterpart of this problem has triggered 
a lot of interest 0, i, la H El in • Ideally, these mves- 
tigations require a well defined system that can be tuned 
all the way from the classical to the quantum limit. 

In this perspective, Josephson phase circuits are ap- 
pealing experimental systems. Their phase dynam- 
ics mimics the one of a fictitious particle oscillating 
in a washboard-like potential well, with associated dis- 
crete energy levels. The system anharmonicity is tuned 
through external bias parameters, such as current or 
magnetic flux. In a two-level system with a sufficiently 
long coherence time, Rabi oscillations appear during the 
transient following the switching on of the microwave 
driving field [13]. A multilevel system also supports sim- 
ilar (Rabi-like) oscillations In this regime, the en- 
ergy granularity is still strong as compared to the mean 
energy. Thus the multilevel dynamics was first fully de- 
scribed using a quantum approach. Later, Gronbech- 
Jensen et al. have reproduced these oscillations using 
a purely classical theory [14]. Moreover, it was shown 
that multilevel Rabi-like oscillations were less sensitive 
to thermal spoiling than their 2-level counterpart [15]. 
These results raise the intriguing question addressed in 
the present paper: are Rabi-like oscillations an unam- 
biguous proof of quantum behaviour? Besides this funda- 
mental motivation, a fine understanding of the multilevel 
regime is also of interest for the accurate control of phase 
qubits, and has motivated recent theoretical [H, U?] and 
experimental works p^, [191, |20| . 

Very recently, it was proposed to analyse the reso- 



nance lineshape under strong driving to track classical 
and quantum signatures [21[. In this paper, we adopt 
a complementary time domain-approach. We report the 
measurement of Rabi-like oscillations in a dc SQUID cir- 
cuit for various /iw amplitude. Improved coherence times 
with respect to previous experiments [1] enable us to ex- 
plore the crossover from two-level to multilevel dynam- 
ics in detail. When up to 3 levels are involved in the 
oscillation, the experimental Rabi frequency versus the 
excitation amplitude exhibits a clear quantum behavior, 
which can definitely not be explained by a classical ap- 
proach. However, when more than 4 levels participate 
in the oscillation, the classical and quantum predictions 
converge. We discuss this issue in the light of Wigner 
function calculations, showing in particular that when 
the system energy reaches its maximum value, quantum 
interference effects are strong, even in the case of a high 
power excitation. 



II. THE DC SQUID AND ITS OPERATION 

The anharmonic oscillator under consideration is real- 
ized in a current biased dc SQUID, which consists of two 
identical Josephson junctions (JJ), each with a critical 
current /q and a capacitance Cq. The junctions are em- 
bedded in a superconducting loop of total inductance L^, 
threaded by a flux <l>5. As discussed in Ref. [22], the phase 
dynamics of the SQUID can be treated as that of a fic- 
titious particle having a mass m = 2Co(<^>o/27r)^ moving 
in a one-dimensional cubic potential (<l>o = h/2e is the 
superconducting flux quantum). The potential is com- 
pletely characterized by the frequency Up of the bottom 
of the well and a barrier height AU [Fig. [1]. These two 
quantities depend on the magnetic flux and vanish at the 
squid's critical current /c. For bias currents 1^ < Ic, 
the particle is trapped in the anharmonic potential well. 
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FIG. 1: Principle of the SQUID operation, (a) Schematic of 
the electrical circuit, (b) Digital sampling oscilloscope record 
of the high frequency flux signal applied on the superconduct- 
ing loop: the /iw excitation pulse is followed by a measuring 
dc pulse which adiabatically reduces the barrier height and 
induces a selective tunnel escape of excited states. 



In the presence of an applied driving magnetic flux, the 
corresponding Hamiltonian reads: 
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Here, oOp = 27TUp] x is the phase along the escape direc- 
tion, and p is the conjugate momentum associated with 
the charge on the junction capacitance, Q = {2e/h)p. 
The cubic coefficient is given by: 
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The /iw flux acts as an oscillating force with a strength 
f^yj and a frequency u. Starting from Eq. ([T]), the SQUID 
dynamics can be treated either classically or with quan- 
tum mechanics. In a quantum description and in the 
absence of driving flux signal, the quantized vibration 
states within the potential well are denoted |n), corre- 
sponding to energy levels En with n = 0, 1, ... The fre- 
quency associated to the |n) \k) transition is denoted 
Unk' The parameter S = uqi — vi2 fully characterizes the 
system anharmonicity. At second order in a perturbative 
approach: 
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is adiabatically switched on through the SQUID at flxed 
magnetic flux $5. The working point (/fo, <l>5) deflnes the 
geometry of the potential well and the circuit is initially 
in the ground state. A fast composite flux signal, pre- 
sented in Fig. [T][b), is then applied. The /iw excitation 
pulse is followed by a dc flux pulse which brings the sys- 
tem to the measuring point within a few nanoseconds. 
This flux pulse reduces adiabatically the barrier height 
AU and allows the escape of the localised states to flnite 
voltage states. Adjusting precisely the amplitude and du- 
ration of the dc flux pulse, it is theoretically possible to 
induce a selective escape of excited states. Because the 
SQUID is hysteretic, the zero and flnite voltage states 
are stable and the result of the measurement can be read 
out by monitoring the voltage Vs across the dc SQUID. 
Our measurement procedure destroys the superconduct- 
ing state and has to be switched off to reset the circuit. 

Though adiabatic, our measurement is one of the 
fastest implemented in Josephson qubits. The delay be- 
tween the end of fiw and the measurement (top of the dc 
pulse) is less than 3 ns. The measurement time, corre- 
sponding to the top of the pulse, is 4 ns. The measure- 
ment speed is crucial to efficiently detect fast relaxing 
states. This is especially true for a multilevel system, 
since the relaxation rate of the excited state \n) scales as 
n. Detection in the multilevel regime has been charac- 
terized in Ref. [23] for settings identical to those used in 
this article. For states described by the set of occupan- 
cies {po,Pi,P2,P3}, the escape probability Pg out of the 
potential well reads: Pe = 0.03 + 0.54pi + 0.79p2 + 0.90p3. 

The SQUID studied in this article consists in two large 
JJs of 15 /im^ area (/q = 1.242 fiA and Co = 0.56 pF) en- 
closing a 350 /im^-area superconducting loop. The two 
SQUID branches of inductances Li and L2 contribute to 
the total loop inductance Lg = 280 pH with the asym- 
metry parameter 77 = (Li — L2)/Ls = 0.414. As in our 
previous work, these parameters were determined com- 
bining two independent sets of measurements: i) macro- 
scopic quantum tunnelling of the ground state and ii) 
spectroscopic study of the |0) |1) transition [23]. 

The chip is cooled down to 30 mK in a dilution fridge. 
The thermal energy is then small compared to the oscil- 
lation energy of the circuit hup ^ x 500 mK. The high 
frequency flux signal is guided by 50 Q coax lines and at- 
tenuated at low temperature before reaching the SQUID 
through a mutual inductance. The nominal room tem- 
perature /iw amplitude is denoted in the following V^w 
For details on the circuit fabrication and the experimen- 
tal setup description, we refer to Ref. [23i] . 



the detuning between the transition |n — 1) \n) and 
the fundamental one |0) |1) scales as z^oi — J^n-i,n = 
S X [n — 1). In the following, for improved precision in 
the data analysis, the eigenenergies En of the system were 
determined numerically. 

Our procedure to perform experiments consists in the 
repetition of an elementary sequence, which is decom- 
posed into four successive steps. First, a bias current 



III. RABI-LIKE OSCILLATIONS 

In this section, we present experimental Rabi-like os- 
cillations and their classical and quantum description. 
Particularly, we focus on the dependence of the Rabi- 
like frequency on the microwave excitation amplitude: 
we examine when this particular measurement offers a 
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ment pulse settings, remain unchanged. During the tran- 
sient regime, Rabi-hke oscihations (RLO) are observed, 
as shown in Fig.[2^c). From a fit of the first oscillations to 
a damped sine function, the RLO frequency ur and the 
damping characteristic time T2* are extracted. It reaches 
here 25 ns, confirming the coherence improvement in this 
sample. The dependence of ur with the jaw amplitude 
is plotted in Fig. [3l The error bar on is smaller than 
the experimental point size, around 5 MHz. In our ex- 
periments, we have chosen to excite levels which do not 
leak to voltage states before the measurement pulse. This 
implies a limitation on the maximum microwave excita- 
tion amplitude. Experimentally, this was checked by per- 
forming an experiment without measurement pulse, and 
by comparing the escape probability with and without 
microwaves. 



B. Classical interpretation 



FIG. 2: (a): low power spectroscopy (VJrxio = 17.8 mV). The 
fit to a Gaussian lineshape gives the resonance position and 
the linewidth. (b): energy relaxation as a function of the 
measurement delay Tm, fitted to an exponential law with the 
damping time Ti. (c): Rabi-like oscillation observed in the 
transient regime for high power /xw driving {V^w =631 mV). 
The oscillation frequency and the damping time T2 are ex- 
tracted from a fit to an exponentially damped sine function. 



clear signature of quantum behavior. The observability 
of beatings, predicted by the quantum approach, is also 
discussed in the end of the section. 



A. Experimental procedure 

Hereafter we present results obtained at the working 
point {h = 2.222 /iA, <I)6 = -0.117 ^0). The experi- 
mental low power spectroscopy measurements (Fig.[2fa)) 
present a resonant peak centered at 8.283 GHz with a 
full width at half maximum Au = 110 MHz. The res- 
onance is interpreted in a quantum description as the 
transition from |0) to |1) at the frequency uqi- Having 
the SQUID electrical parameters in hand, the character- 
istics of the potential well are calculated. The plasma 
frequency is Up = 8.428 GHz and the well contains 8 en- 
ergy levels for a total depth AU = h x 61.6 GHz. Com- 
pared to Ref. [l|, the present sample has lower decoher- 
ence {Au = llOMHz instead of 180MHz), and the chosen 
bias point displays higher anharmonicity {6 — 160 MHz 
instead of 100 MHz). This improvement allows for much 
more transition-selective coherent /iw driving. 

The transient non linear dynamics of the SQUID is 
probed using a fast /iw flux pulse followed by a measur- 
ing dc pulse. The /iw excitation is tuned on the reso- 
nance frequency obtained from low power spectroscopy 
and the /iw pulse duration is increased from 2 ns to 80 ns. 
The measurement delay, as well as the other measure- 



To describe the observed Rabi-like oscillations, we first 
apply the classical approach developed in Ref. [14]. Hav- 
ing turned on the resonant driving force, the particle en- 
ergy oscillates during a transient regime. The underly- 
ing mechanism is related to the potential non-linearity, 
which induces a dephasing between the particle oscilla- 
tions and the /iw excitation. It makes the /iw alterna- 
tively accelerate and brake the particle oscillations within 
the well. In the general case, vr depends on f^^^, v and 
T\ , and the complete set of equations in Ref. ^J] must be 
solved. Here, the excitation was tuned on the resonance 
frequency obtained from low power spectroscopy, which 
corresponds to the bottom well frequency in a classi- 
cal approach. In the limit of small damping (vrTx ^ 1) 
which corresponds to the experimental situation, vr is 
then insensitive to T\ and coincides with the value cal- 
culated in the absence of dissipation. In the following 
analysis, we solved numerically the equation in Ref. [ij] 
under these two conditions. For the particular situation 
of small iiw excitation {f^w/ci <C 1), an analytical ex- 
pression can be derived: 



4 ^ ^ \muL)^J 

where A is a dimensionless parameter given by: 
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Fig. [31(a) shows the fit of experimental data to the clas- 
sical theory, the calibration coefficient c = f^w/V^w be- 
ing the only free parameter of the model. For consis- 
tency, we should note that, in a classical description, the 
low power resonance presented in Fig. [2ja) is associated 
with the plasma frequency characterizing the bottom of 
the well. As spectroscopic data are used to extract the 
SQUID parameters, the classical analysis leads to slightly 
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FIG. 3: (Color online) Fit of the experimental vr (•) to a classical (a) and quantum theory (d). In the case of weakly damped 
Rabi oscillations, vr is independent of the damping time and identical to the one of a lossless system, both in the quantum 
and classical cases. The classical model (solid line in (a)) fails to describe the low amplitude regime, unless excessive losses, 
incompatible with the measured Ti are introduced (dashed line). The quantum prediction (solid line in (c)) accurately describes 
the experiment over the whole excitation amplitude range. For vi > z^oi — z^02, the curve deviates from the two level theory 
(dashed line), indicating a multilevel dynamics. Simulated Rabi oscillation for a /xw amplitude corresponding to Fig. [2jc) are 
also shown under various hypothesis: classical calculation with Ti = 60 ns (b) and Ti = 2.5 ns (c) and quantum prediction 
without (e) and with decoherence (f) which suppresses beating. 



different electrical parameters; they are listed in note [2^. 
In the weakly damped limit, the best adjustment which 
is obtained for c = 3.50 x 10~^^ J.V~^ does not describe 
correctly the data on the full range of excitation ampli- 
tudes. At low fiw amplitude, the observed linear depen- 
dence is strongly different from the /^^^ law predicted 
by the classical approach. However, the agreement im- 
proves for higher driving amplitude; this convergence will 
be discussed at the end of the paper. 

As a remark, it is possible to fit the data if arbi- 
trary high losses are introduced. Decreasing Ti down 
to 2.5 ns and using c = 2.90 x 10~^^ J.V~^, one obtains a 
good agreement between the data and the classical the- 
ory (dashed line in Fig. [3]^a)). However, this hypothesis 
is inconsistent with the measured Ti = 60 ns. It also 
would lead to overdamped Rabi-like oscillations, with a 
T2 = 3 ns (Fig. [31(c)) which is one order of magnitude 
smaller than the measured one. 



C. Quantum approach 

1. Decoherence- free calculation 

The explanation of experimental data requires a full 
quantum treatment. Since moderate decoherence has a 
minor impact on the Rabi-like frequency, we first develop 
a pure Hamiltonian theory [l|, [Tol, [27t . 



The particle is assumed to lay initially in the ground 
state. At t = 0, the /iw flux is instantaneously switched 
on, and the particle quantum state |^(t)) evolves into 
a superposition of the level states |n) according to the 
Schrodinger equation with the Hamiltonian: 



H{t) hup 
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(6) 



where x and p are the reduced phase and momentum 
operators defined by: 
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The expression Eq. (|6|) contains all the characteristic fre- 
quencies of the problem: i) the bottom well frequency 
Up^ ii) the detuning 5 = vqi — V12 linked to the potential 
anharmonicity and given by Eq. (|3]) and iii) z^i, the Rabi 
frequency in the 2-level limit for a microwave excitation 
frequency v = i^oi, which reads: 
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To solve this time-dependent problem, we introduce a 
new picture in which the evolution of the particle state 
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I^Q(t)) is driven by the Hamiltonian HQ^t). This new 
picture is related to the Schrodinger one by: 



{n|*oW) 



(9) 



In the 2-level hmit, this transformation corresponds to 
the introduction of the weh-known rotating frame at fre- 
quency u. To proceed, we apply the rotating wave ap- 
proximation (RWA) to simplify the Hamiltonian HQ(t). 
Considering a potential well with trapped energy lev- 
els, the resulting Hamiltonian Hq^rwa can be expressed 
in the {\n)} basis as follows: 



0,RWA 



= h 



( t 



V 











(10) 



Here, lS.n(y) — i^on — and the ground state energy Eq 
has been set equal to 0. Note that x presents matrix el- 
ements of the order of y^S/up on the diagonal and other 
which connect neighbor states of order 2. Other terms 
of higher order in \/Sji^ are also present in the matrix. 
These terms vanish in the RWA, which remains valid as 
long as ui <C i^oi. In this quantum approach, the experi- 
mental low power resonance is attributed to the |0) |1) 
transition; thus the following analysis is performed for 
u = uqi. The eigenvalues associated to the eigenstates 
\eno) of ^o,RWA are shown in Fig. [H Having solved the 
time evolution of starting from the initial state 

|0o) = |0), the conversion to the Schrodinger picture is 
done according to Eq. (|9]). Particularly, the probability 
to find the system in the state |n) in the course of the 
oscillation reads for n > 0: 

N-l 

Pn{t) = 2 ^ Kn|e/eo)(e/co|0)(0|e^o)(e/ol^)l 

k,l>k ^^^^ 

X cos[27r(Afe — Xi)t]. 

The occupation dynamics of state |n) implies the fre- 
quency set {|A/c — A^|} (/c 7^ /) with weights associated to 
the connection of |0) to |n) via the eigenstates of Hq^rwa- 
As shown in Fig. [3l^e), the average energy of the parti- 
cle (H) = '^nPn{t)En Undergoes Rabi-like oscillations. 
The model also predicts modulations of the RLO am- 
plitude which will be discussed in the next paragraph. 
Our analysis shows that an estimate of ur can be ob- 
tained by taking the minimum value of the frequency 
set {\Xk — A/|}. For improved precision, the theoretical 
Rabi-like frequency is extracted from a fit of the first os- 
cillations to an exponentially damped sine function. This 
approach was checked to give similar results as a direct 
numerical integration of the time dependent Schrodinger 
equation. 

As shown in Fig. [3l[d), the quantum model describes 
all the experimental features of the versus f^^ depen- 
dence with c = 3.05 X 10~^^ J.V~^. The set of experi- 
mental data covers the crossover between the 2-level and 
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FIG. 4: (Color online) First eigenvalues of Hq^bwa as a func- 
tion of ui. The dashed lines correspond to the limit 1^1/ 6 <^ 1. 
The inversion of the spectrum is due to the rotating frame 
transform. 



the multilevel dynamics. Indeed, when the /iw ampli- 
tude is small compared to the anharmonicity {ui <C ^), 
we have |eoo) = (|0) + \l))/V2 , |eio) = (|0) - \1))/V2 
and \eno) = 1^) foi" the other higher energy levels. The 
dynamics thus concerns the first two levels and one re- 
trieves the familiar result ur = vi [see Fig. |4]. This is 
nearly the case for the first measured point, which corre- 
sponds to vi = 65 MHz and presents a contamination of 
level 1 2) below 10%. When ui ^ the coupling between 
neighboring levels starts to distort the energy spectrum 
of the lowest energy levels of Hq^rwa [Fig. [4]. Conse- 
quently, the state |e2o) gets contaminated by |0) and |1) 
and the dynamics involves now 3 levels. In this regime, 
the ur dependence on ui starts to deviate from a linear 
behavior, indicating the onset of two-photon processes 
[19, 20]. When ui > an increasing number of levels are 
involved. As an example, the experimental oscillation 
presented in Fig. [2](c) involves 4 states. At larger /iw 
amplitude, a macroscopic number of levels are involved; 
this regime is not achieved at this working point but was 
discussed in Ref. [1] with about 10 levels involved. The 
multilevel dynamics is characterized by a clear saturation 
of the ur dependence, compared to the low power linear 
behavior. 



2. Beating and decoherence 

In the absence of dissipation, the quantum model pre- 
dicts a modulation of the multilevel RLO amplitude. We 
interpret this phenomena as a reminiscence of the dif- 
ferent frequencies involved in the system dynamics. The 
modulation becomes more and more complex with in- 
creasing /iw amplitude. Indeed, the higher the ampli- 
tude, the larger the number of involved levels and the 
larger number of frequencies in the system. A clear beat- 
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ing appears when the frequency set {|A/e — A/|} displays 
two minimal values close to each other. 

However, these modulations were not observed in our 
experiments. We first analyzed the possibility of beat- 
ing suppression caused by the finite rise time of the /iw 
amplitude. If it is too slow, the population of the states 
le^o) for n > 2 can be severely decreased as compared to 
the instantaneous hypothesis. In the limiting case of an 
adiabatic increase of the /iw amplitude, the particle state 
remains in the {leoo)? kio)} subspace. |Ai — Ao| is then 
the only frequency involved in the system dynamics (see 
Fig. Hj) and there is no beating. For a rise time as large 
as 2 ns (twice the experimental rise time), we checked 
numerically that this effect does not lead to destructive 
interference strong enough to suppress the beating. 

We therefore argue that the strong suppression of the 
beating phenomenon could be due to decoherence. The 
dissipative dynamics of the oscillator was investigated us- 
ing a phenomenological model. We assume that the re- 
duced density matrix of the driven oscillator p satisfies 
the master equation 

ihdtp = [^o,RWA, p] - ijCnlp] - ^^d[p] - (12) 

Here, Cji describes relaxation (accompanied by dephas- 
ing) and can be expressed as follows: 

^1z[p\kn = ^ [(7/cm + 7nm) Pkn " ^^^knlrnkPrnm] • (13) 
m 

The coefficients jkn can be calculated microscopically in 
the Markovian limit, assuming the environment is a bath 
of harmonic oscillators. The last term Cx) describes the 
pure dephasing: 

^v[p]kn = (1 - Skn)^knPkn , (14) 

with Akn = ^nk ^ 0. This term is phenomenological: 
the coefficients A were obtained microscopically only in 
the two-level limit. According to our analysis, relaxation 
alone (corresponding to Ti = 60 ns in the experiment) 
is too weak to explain a strong suppression of the beat- 
ing. Instead, the absence of beating can be attributed to 
strong pure dephasing between next-neighboring levels, 
described by An,n+2- Indeed, the observed exponential 
decay of the Rabi-like oscillations with the suppressed 
beating can be well explained assuming A02 ^ ha/vi^ 
which becomes relevant for the measured escape proba- 
bilities corresponding to the cases of vi = 130 and 260 
MHz with a = 5 X 10^ GHz^. 



IV. DISCUSSION 

Surprisingly, when vi > 1.5(5, the classical and quan- 
tum predictions for the i'r{i'i) dependence are quite iden- 
tical. For the lowest excitation amplitudes meeting this 
condition, only 4 levels are involved in the Rabi-like oscil- 
lation: the granularity of energy as compared to the mean 



energy of the system is clearly not negligible. Moreover, 
this convergence is obtained for weakly damped oscilla- 
tion, for which moderate relaxation and dephasing do 
not affect vr. Thus, this effect is not due to decoherence 
which would degrade the system quantum dynamics into 
a classical one. The convergence is linked to the intrinsic 
physics of a weakly anharmonic oscillator. 

To shed some light on the interplay between non- 
linearity and the quantum or classical nature of the Rabi- 
like oscillations, we developed an approach similar to the 
one presented in Ref. [ll|. We calculated the time evo- 
lution of the Wigner function 

W{x,p,t) = \ I dx'{x^x'\p{t)\x-x')e-^'P'''^^, (15) 

Tih J 

during one Rabi oscillation, starting from the oscillator 
ground state |0). Here, p is the reduced density matrix 
describing the driven oscillator. For the sake of compar- 
ison, we also determined the classical time evolution of 
the probability density for a gaussian distribution of ini- 
tial points in phase space which mimics the distribution 
associated to the ground state |0). Decoherence effects 
were neglected in both approaches. This analysis high- 
lights the differences in the system dynamics arising from 
quantum mechanics. Fig. [5] shows the results for three 
different excitation amplitudes: i) the lowest excitation 
amplitude measured in this paper, vi = 60 MHz which 
induces a quasi 2-level oscillation, ii) vi = 260 MHz, the 
highest excitation amplitude and iii) vi = 520 MHz, a 
(non measured) value which generates 7-level oscillation, 
similar to the high excitation results presented in Ref. . 

The red areas in Fig. [5] correspond to negative val- 
ues of the Wigner function, an unambiguous signature 
of the quantum nature of the state under consideration; 
in the blue regions W is positive. As expected for short 
times (shorter than the so-called Ehrenfest time), the 
system behaves classically and the positive part of W re- 
sembles the classical probability distribution, regardless 
of the value of the microwave amplitude, as illustrated 
in Fig. [5] for t = O.ITr. For longer times, this resem- 
blance gets lost for small microwave amplitudes. This 
is to be expected: due to the anharmonicity, the system 
remains close to the two-level limit for small amplitudes 
and the time evolution is essentially quantum. For larger 
amplitudes, the dynamics involves an increasing number 
of quantum states and according to the correspondence 
principle should behave more and more classically. In- 
deed, the time evolution of the blue regions closely fol- 
lows that of the classical probability distribution in this 
limit. 

Interestingly, though, is the fact that the red regions of 
negative W continue to exist at any value of microwave 
amplitude. These regions, signalling the occurrence of 
quantum interferences, appear to be most pronounced 
for t = 0.5Tr. a qualitative explanation for this fact can 
be obtained by considering the time evolution of the an- 
harmonic oscillator, starting from its ground state |0) at 
time t = 0. Under the influence of the microwaves, reso- 
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FIG. 5: Quantum (q) and classical (c) dynamics during the first Rabi period Tr, for increasing excitation amplitudes ui. 
Decoherence is neglected in all calculations. Wigner and classical distributions are plotted in the phase space (position x, 
conjugate moment p). In the Wigner distribution plots, red (blue) corresponds to negative (positive) values. The ticks indicate 
the position of the origin of the phase space. Negative areas are always intense close to t = 0.5Tr, when the system experiences 
fully the non linearity of the potential well, highlighting the quantum nature of the corresponding states. Note also the 
progressive convergence, as vi increases, between the principal blue cloud of the Wigner map and the classical prediction. 



nant with the transition between the lowest eigenstates, 
transitions between neighboring levels are possible, and 
a superposition state is formed involving an increasing 
number of energy levels. Due to anharmonicity, the de- 
tuning increases with increasing level index and the tran- 
sitions become less resonant. At t = 0.5T/^, non-linearity 
prevents higher energy states from being excited and mi- 
crowaves lead to a reduction of the system energy. Thus 
the existence of Rabi-like oscillation is closely linked to 
the anharmonicity of the potential well and the max- 



imum number of states involved in the superposition is 
determined by the interplay of the anharmonicity and the 
microwave- amplitude. Such oscillations do not show up 
in a harmonic oscillator driven by resonant monochro- 
matic excitation. In this last situation, the generated 
quasiclassical states present minimum fluctuations and 
are thus very close to classic ones. On the contrary, in 
the case of an anharmonic oscillator, for times close to 
t = 0.5Ti^, non-linearity gives birth to pronounced quan- 
tum interferences, and the corresponding Wigner func- 
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tion is characterized by the existence of large negative 
regions. 

Finahy, we can remark that the analysis of the Rabi 
oscillation frequency is one measurement among many 
possibilities. It appears that this quantity is not very 
sensitive to the quantum interferences characterising the 
extreme state of the Rabi oscillation. To track a quan- 
tum signature in the multilevel regime {ui > J), one could 
concentrate on the oscillation beating, a reminiscence of 
the different transition frequencies in the system. An- 
other appealing perspective would be to extend the state 
tomography scheme developed in phase Qubits [26] to 
our multilevel system. Though experimentally challeng- 
ing, this would open interesting possibilities such as the 
investigation of the effects of decoherence on the state 
superposition generated at high excitation power. 

V. CONCLUSION 

In conclusion, a dc SQUID is a tunable and well de- 
fined anharmonic quantum oscillator which appears as 
an experimental model system to explore quantum and 



classical dynamics. The dependence of the Rabi-like os- 
cillation frequency on the jaw amplitude is a clear probe 
of the nature of oscillator dynamics when up to 3 levels 
participates in the oscillation. It then exhibits a clear 
quantum behavior which is not described by a classical 
model. At higher excitation power, when more than 4 
states are involved in the oscillation, the quantum and 
classical predictions for this particular quantity converge. 
This convergence was discussed with the help of Wigner 
quasi-distribution of probability. However, one should 
not hastily conclude than this non linear oscillator be- 
haves classically as soon as 4 levels are excited. Indeed, 
our theoretical study shows that for times corresponding 
to the middle of the Rabi-like oscillation (the general- 
ization of a TT pulse), the system always presents strong 
quantum interferences, even in the high power excitation 
regime. Future work will concentrate on the shape of 
the oscillations which contains more information on the 
multilevel dynamics. 

We thank F. Faure, W. Guichard, L. P. Levy, and A. 
Ratchov for fruitfull discussions. This work was sup- 
ported by two ACI programs, by lUF and IPMC and 
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